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Abstract: We explore the high and low temperature behavior of non-local observ- 
ables in strongly coupled gauge theories that are dual to AdS. We develop a systematic 
expansion for equal time two-point correlation, spatial Wilson loops and entanglement 
entropy at finite temperature using the AdS/CFT correspondence, leading to analytic 
expressions for these observables at high and low temperature limits. This approach 
enables the identification of the contributions of different regions of the bulk geometry 
to these gauge theory observables. 
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1 Introduction 



The finite temperature behavior of non-abelian gauge theories has been of interests 
to physicists over the years (see [1-11] and references therein). Of particular interest 
is the behavior at strong couphng where we have a powerful tool in the AdS/CFT 
correspondence [12-15]. This correspondence has been used in various areas ranging from 
Relativistic Heavy Ion Collider (RHIC) to condensed matter physics. The study of non- 
local observables (e.g. two-point function, Wilson loop, entanglement entropy) at finite 
temperature for strongly coupled gauge theories using the AdS/CFT correspondence is 
an interesting problem in its own right. One of the questions we focused on, is the regions 
in the bulk that contribute most significantly to these observables of the 'boundary' 
gauge theory in the different limits. This led us to construct analytic expansions both 
at high and low temperature. 

In this paper, we will consider strongly coupled large-N gauge theories in d-dimensions 
that are dual to AdS^+i. As mentioned above, our goal is to understand analytically the 
behavior of non-local observables like equal time two-point function, spatial Wilson loop 
and entanglement entropy in the low and high temperature limits. We will henceforth 
develop a systematic expansion using the AdS/CFT dictionary. All of these observables 
of the boundary gauge theory, associated to a region A can be computed by calculating 
the area of the bulk extremal surface anchored on the boundary of A (for the two-point 
function the relevant quantity is the geodesic, for what follows we will use the label 
'extremal surface' for the sake of simplicity). As a consequence of conformal invariance 
the dimensionless parameter of interest is Tl, when we are computing some observable 
with length scale / at temperature T. 

We develop separate expansion techniques at low and high temperature. At low 
temperature, i.e. T <C l/l, the extremal surface is restricted to be near the boundary 
region and thus the leading contribution to the area comes from the AdS-boundary. This 
contribution is just the zero temperature result which has been thoroughly investigated 
and well understood in the literature. Finite temperature corrections correspond to the 
deviation of the bulk geometry from pure AdS. At low temperature, these corrections 
are small and can be computed perturbatively. On the other hand, at high temperature, 
i.e. T ^ l/l, the physics is completely different and more interesting. As the thermal 
fluctuations become more and more significant, the extremal surface associated with the 
observable approaches the horizon.^ At high temperature, the extremal surface tends 
to wrap a part of the horizon and the leading contribution comes from the near horizon 
region of the surface. This has following consequences: 1) the equal time two-point 
function decays exponentially; 2) the leading contribution to the entanglement entropy 

"'^As recently discussed in [20], extremal surfaces of any dimensionality and anchored on arbitrary 
shaped region on the boundary at constant time can not penetrate the horizon in any static spherically 
symmetric space times. At high temperature, we will show that the extremal surface approaches the 
horizon exponentially fast but it always stays finite distance above the horizon. 
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(spatial Wilson loop) is proportional to the volume (area) of the associated region. The 
sub-leading terms are more complicated because they receive contributions from full bulk 
geometry. The sub-leading terms in various observables, in particular the entanglement 
entropy, may possibly contain important information about these gauge theories. 

Fascinating question arises when one considers the high temperature regime of the 
field theory. At weak coupling, it is known that the high temperature behavior of the 
gauge invariant correlation functions (e.g. (trF^y(Z)trF^^(O)), spatial Wilson loop) is 
given by classical statistical mechanics[l], this obviously does not apply to thermody- 
namic quantities like the specific heat. Indeed, one expects on the gauge theory side, 
that the leading contributions to the gauge invariant correlation functions with depen- 
dence on one length scale / come from modes with wavelength ~ /. As the system is 
heated the modes with wavelength ~ / become more and more populated. Therefore, 
equal time two-point functions and the spatial Wilson loops should behave classically 
at high temperature (T ^ 1//). It is tempting to speculate that at strong coupling an 
analogous statement can be made. This raises a puzzling question. On the one hand, 
the high temperature behavior is given by the near horizon geometry and on the other 
hand we have speculated that it is given by classical statistical mechanics. Can we then 
conclude that the near horizon contribution corresponds to classical statistical mechanics 
on the boundary? 

In contrast, we found that it is easier to understand the high temperature behavior 
of entanglement entropy, which is a fully quantum mechanical concept and does not have 
a classical analog. At high temperature (T ^ 1//), the leading finite contribution to 
the entanglement entropy comes from the near horizon region and it is just the thermal 
entropy. Whereas, contributions arising from deviations away from the near horizon 
region are responsible for the sub-leading terms, which actually measure quantum en- 
tanglement between the region and its surroundings. Possibly this is just a manifestation 
of UV/IR duality [16]. 

The rest of the paper is organized as follows. We start with a brief review of how the 
temperature of the boundary gauge theory is related to the horizon of the bulk theory in 
section 2. In section 3, we develop a systematic expansion for the equal time two-point 
function at finite temperature. Then using that expansion we examine the high and 
low temperature behavior of the two-point function. In section 4 and 5, we go on to 
develop similar systematic expansions for rectangular spatial Wilson loop and entangle- 
ment entropy of a rectangular strip respectively at finite temperature. Consequently, we 
obtain analytic expressions for the spatial Wilson loop and the entanglement entropy 
at high and low temperature limits. Finally, we summarize all the results in section 6 
by considering the prototype case of four-dimensional A/" = 4 Super Yang-Mills gauge 
theory. Then we conclude in section 7 with future directions. Several technical details 
have been relegated to two appendices. 
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2 Gauge theories at finite temperature: Schwarzschild-AdSd+i 



In this paper, we will consider strongly coupled large-N gauge theories in c?— dimensions 
that are dual to AdS^+i. At nonzero temperature, all the thermal effects can be analyzed 
simply by introducing a black hole in the bulk. In other words, at finite temperature 
AdSrf+i should be replaced by Schwarzschild-AdS^+i 



where, R is the AdS radius. The temperature T is given by the Hawking temperature of 
the black hole which can be determined by demanding that the Euclidean continuation 
of the metric (2.1) 

^2 / j,d \ ^2 p2 

ds' = -Jl-^]dtl + -^dx' + — ^dr^ (2.2) 



r2 ( 1 - ^ 



is regular at the horizon. Near the horizon r = rn, the Euclidean metric looks like 



^2 

C^Snear horizon = p'^d(j)^ + dp^ + (2-3) 



where p and are defined as 



R^jr-rn) rnd 
^ = 'V r^d ' '^=2R-'^- ^'-'^ 

Now (f) must be periodic with period 27r in order to avoid a conical singularity at p = 
(i.e. r = th)- From equation (2.4) it is clear that Euclidean time t^; is also periodic and 
the period can be identified with 1/T. Therefore, the temperature of the (i-dimensional 
boundary field theory (after restoring K) is given by, 

^ h rud 

3 Two-point function 

We will compute the equal time two-point function of some gauge-invariant scalar op- 
erator 0{t,x) in the limit ^ 1, A = Oym^ ^ 1- According to the AdS/CFT corre- 
spondence, in this limit, the calculation of two-point function [13, 14] of gauge-invariant 
operator 0{t,x) simply reduces to solving classical equation of motion for some bulk 
field that couples to 0{t, x). This prescription has been successfully generalized [17, 18] 
for the finite temperature case (for a nice review see [19]). 

For our purpose it is more convenient to use the alternate prescription proposed in 
[21]. Following [22] the equal-time two-point function {0(t, x)0(t, y)) can be represented 
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as a path integral that sums over all the paths that join the boundary points {t, x) and 

{0{t,x)0{t,y)) = J DPe-^^(^), (3.1) 

where L{V) is the proper length of the path and A is the conformal dimension of the 
operator O. For operators with large conformal dimension, we can perform a saddle 
point approximation^ 

{0{t,x)0{t,y))^ ^ e-^^ (3.2) 

geodesies 

where £ is the geodesic length of the geodesic between boundary points {t, x) and (t, y). 
The geodesic length C has a divergence that comes from the boundary. This divergence 
can be removed by introducing a cutoff ri, and then defining renormalized geodesic length 
by 

£ren = - 21nr;,. (3.3) 

Note that the divergent piece is independent of both temperature and dimension. Finally, 
the renormalized two point function is given by, 

(C(t,a;)C(t,?/)) ^e"^^-". (3.4) 

3.1 Finite temperature expansion 

The bulk-metric is given by'^ 

2 / 1 ' H 



ds^ = -r^ ( 1 - -f ) dt' + r^dx^ + — ^dr\ (3.5) 

( 1 — - 



We can choose our coordinates so that the two points are (t,x = — 1,0, ...) and (t, x 



2' 

|, 0, ...) respectively. Therefore, the relevant part of the Schwarzschild-AdSrf+i metric is: 



\ 

ds^ = r'^dx^ + — ^dr\ (3.6) 

^2 M _ ILff 



In the affine parametrization, the geodesic equations are given by. 



r 



x=-^, (3.7) 



^2 ^ tytd 



j \ lyd 

^It is important to note that for Lorentzian correlators this approximation should be used more 
carefully [23-25]. 

^For the sake of cleanliness, we will use AdS radius R= \. We will restore R by dimensional analysis 
whenever necessary. 
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where, is an integral of motion associated with the Kilhng vector and r = 
represents the point of closest approach of the geodesic. We also have taken the affine 
parameter to represent the geodesic proper length s. It is clear from equation (3.8) 
that each geodesic has two branches, which we will denote as a;+(r) and x_(r), joined 
smoothly at (r = Tc^x = 0). In principle, rc can be determined using the boundary 
conditions: 

x_(oo) = --, x+(oo) = -. (3.9) 
Now using equations (3.7, 3.8), for the positive branch of the geodesic, we obtain 

dr r'^ I f r1\ f r'ir . , , 




Integration of the last equation leads to 



/ /"°° r^dr I ^ r'jj 



d \ -1/2 



^d 



udu 



_H_ 

Tc Jo V 

Unfortunately this integration can be performed analytically only for Schwarzschild- 
AdSs (i.e. d = 2). For a general d, we can do a systematic expansion. 

_ 1 f. r[i + n]r[i + f] (r^y^ 

2r.^ r[l + n]r[i + f] \rj ' 

We should be more careful about the convergence of the series (3.13) before we use it. 
For large n, the series goes as ~ ^ {th/^cY'^ and hence the series converges for th/tc < 1. 
For any finite temperature, it can be shown [20] that Vc > tr- Therefore, the sum (3.13) 
is well-defined. 

Next we will calculate the regularized geodesic length by using equation (3.8) 

r dr ( rfjY^^^ 



1/2 

u"" ) . (3.11) 



Where, the factor of 2 comes because of the two branches of the geodesic. The equa- 
tion (3.14) has a divergence that comes from the upper limit of the integration. This 
divergence can be removed by introducing a cutoff and then using equation (3.3) to 
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define the renormalized geodesic length. Now we will proceed to develop a systematic 
expansion for H^en 



Cren=2 J 1 - -^W - 2 In ^ (3.15) 

-2C ^^Vi:Ii±ir!i^V'.--21nn (316) 

=2/' ^^^^ 21nr. + 2r^^f ^ii^f^V'."^ 

(3.17) 

where, we have used the fact that only n = term is divergent. Finally, we have 

i:,en=21n - +> . n wi — • 3.18 

Again it can be shown that for large n, the infinite series in equation (3.18) goes as 
~ n ^^al'^cf'^ ^^"^ hence the series converges for th/tc < 1. In principle, the rest of 
the procedure is very simple; we have to solve equation (3.13) for Tc and then use that 
Tc in equation (3.18) to get the renormalized geodesic length.^ The renormalized two- 
point function is then given by equation (3.4). But in practice, this procedure can be 
performed exactly only for d=2. For 7^ 2, it is not possible to solve equation (3.13) 
analytically to find Tc as a function of /. However at low temperature (i.e. thI ^ 1), we 
can compute Cren perturbatively using equations (3.13, 3.18). And more interestingly, 
equation (3.18) can also be used to determine the high temperature (i.e. thI ^ 1) 
behavior of the two-point function. 



3.2 Two-point function: CFT in (1 + 1) dimensions 

A trivial example, as mentioned earlier, is Schwarzschild-AdS3 (i.e. d = 2). For d = 2, 
equation (3.13) reduces to 

/ 1^ 1 /r^^'" tanh-Wri. 



-) = — (3.19) 

n=0 

Therefore, 



'''' (3.20) 



tanh (^) 



*We would like to stress that the equations (3.13, 3.18) are valid for any temperature. 
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And equation (3.18) becomes 



^rp.rt. 2 In 



E 

n=l 



1 / 



n \ r, 



2n 



=2ln(-1-ln(l 



-In 



H 



2 (cosh(r///) - 1) 



Therefore, 



2A 



{0{t,x)Oit,y)) = rjj 



1 



2 (cosh(rH/) - 1) 
where, / = |x — ?/|. Now using equation (2.5) we finally have 



(3.21) 



(3.22) 



{0{t,x)0{t,y)) = {nT) 



2A 



cosh {2'kT\x — y\) — 1 



3.3 Low temperature two-point function 



(3.23) 



At low temperature, Vc ^ rfj and hence the leading contribution to the geodesic length 
comes from the boundary. The boundary is still AdS and we should get the zero tem- 
perature two-point function as the leading term. Finite temperature corrections can 
be computed by considering deviations from the boundary geometry. At low tempera- 
ture (T/ ^ 1), the corrections to zero temperature result are small and hence can be 
computed perturbatively. From equation (3.13), keeping only few subleading terms, we 
get 



o , Ar(| + 1) frnY , 3v/^r(rf+l 
2r(^) \rj + 8r(rf+|) 



2d 



+ 



(3.24) 



Solving the last equation perturbatively and then using that solution in equation (3.18), 
we get a perturbative expression for Cren (details of the calculation are relegated to 
appendix A). Finally using (2.5), for the rf-dimensional boundary theory at low temper- 
ature we obtain 



{0{t,x)0{t,y)) 



1 



\x-y 



2A 



1+Ci 



2'kT\x — y\ 



d 



+C2 



2'kT\x — y\ 
~d 



2d 



o 



2'kT\x — y\ 
d 



3d" 



(3.25) 
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where 



Ci = - 



¥Ar(|) 



Co 



1 

"64 



vrA {d^ + 2A) 



12v^Ar(rf) 



(3.26) 
(3.27) 



Note that Ci is negative as expected indicating a decrease in two-point correlation be- 
cause of thermal fluctuations. 

Particularly for d = 4, using equations (A.5-A.7), at low temperature limit {T\x — 
y\ <^ 1) we have 



{Oit,x)Oit,y)) 



\x 



|2A 



1 - 



2A fnTlx 



y\ 



15 V 2 
2A(7A + 26) fnTlx 



1575 



7cT\x — y\ 



12 



(3.28) 



3.4 



High temperature two-point function 

At high temperature (i.e. Tl ^ 1), it is more difficult to do a systematic expansion for 
the two-point function. Nevertheless, we can find out the asymptotic behavior of the 
two-point function easily. At high temperature, the leading contribution comes from the 
near horizon part of the geodesic and this piece is easy to compute. On the other hand, 
the full bulk contributes to the subleading terms and hence they are more complicated. 
At very high temperature, Vc approaches r^- The obvious guess of taking the limit 
— ^ rn in equation (3.18) does not work because the infinite series in equation (3.18) 
converges only when Tc < th- But we can rewrite equation (3.18) in a way that allows 
us to take the limit — )■ th without encountering any divergence. 



^ren 2 In 



2\ r[i+n]r[f] frn^^" 



r. . 

n 



^r[i + n]r [l + f] 



2\ ^ fl + nd\T['^+n]T[l + f] /r^" 



n=l 



-^^^ -]+y ^] X V - • (3-29) 



Now, using equation (3.13), we obtain 

£ -2 m V (r / - 2) + f m r + .] r [1 + f ] /^x - 



(3.30) 



The term —2 in the last equation comes from n = term of the series (3.13). The 
infinite series in equation (3.30) converges even for = th and hence the limit — > 
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Boundary 




Horizon 



Figure 1. At high temperature (r/// ^ 1), the actual geodesic (sohd blue line) can be 
approximated by the dashed red line curve that consists of x = — //2,r = th^x = 1/2. 



exists. At high temperature, ~ th and now the leading behavior can be determined 
by taking the limit Tc — )• th in equation (3.30) 



r 

'~'re 



2 In 



th J 



\ , , ^ / 1 ^ r ^1 r [1 + f 1 
+ M-2) + $: - J L .J 



(3.31) 



r[i + n]r[| + f] • 

Therefore, the high temperature two-point function is approximately given by, 

(0(t, x)0{t, y)) ^ r^f e-^'--'. (3.32) 

Where, the prefactor ^^^a is a constant that depends on the dimension d and the con- 
formal dimension A of the operator O 



A, 



■A 



4 



exp 



2-E 



1 \T[l + n]T[l + f] 



^^\ndj r[i + n]r[| + f] 



(3.33) 



Finally, replacing th by corresponding temperature, for the rf— dimensional boundary 
theory we obtain 



{0{t,x)0{t,y)) ^ AdA 



AttT 



2A 



-47rAT\x-y\/d 



(3.34) 



The exponential decay of the two-point function at high temperature can be under- 
stood easily by looking at the geodesic for rnl ^ 1. The actual U-shaped geodesic can 
be approximated by a curve that consists of x = — //2,r = rH,x = 1/2 (see figure 1). 
As thI is increased, the actual geodesic approaches the approximate one. But it can be 
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shown that the proper length of this approximate curve is always greater than the actual 
geodesic and even in the limit r^l — ?■ c>o, two does not coincide. However in the high 
temperature limit, the most dominant contribution to the geodesic length comes from 
the near horizon part which can be reckoned from the approximate curve. The length 
of the near horizon part of the approximate curve is 5 ~ thI- Therefore, it is expected 
that the two-point function {0{t, x)0{t,y)) ~ e~^^"^. 

We have computed the leading behavior of the high temperature {Tl ^ 1) two-point 
function by utilizing the fact that at high temperature, ~ r^- It is a good exercise to 
figure out exactly how close to the horizon the geodesic can reach. As shown in appendix 
B, rc approaches rn exponentially fast 



rn [l + Sd e~ 



(3.35) 



where £d is a numerical constant given by. 



£d 



d 



exp 




r[| + n]r[i + f] v/2 
r[i + n]r[| + f] 




(3.36) 



The geodesic always stays finite distance above the horizon which is consistent with 
[20]. Using equation (3.35), we can calculate the next order correction to the two-point 
function (for details see appendix B). However, the subleading term is exponentially 
suppressed 



{0{t,x)0{t,y))=Ad 



^ (^) 



2A 



-i-KA.T\x-y\/d 



1 + J -A£d e-^^^^'^\^-y\ + 
V d 



(3.37) 



where the dots represent the higher order correction terms. Therefore, at high temper- 
ature, the two-point function is well approximated by equation (3.34). 



4 Spatial Wilson loop 

Wilson loops are another set of important gauge-invariant non-local observables in any 
gauge theory. In a gauge theory, the Wilson loop operator is a path ordered contour 
integral of the gauge field 

WiC) = ^Tr(Ve^c^y (4.1) 

where the trace is over the fundamental representation and C denotes a closed loop in 
spacetime. Expectation values of the Wilson loops are useful to understand the non- 
perturbative behavior of non-Abelian gauge theories and have important applications to 
confinement/deconfinement transitions and quark screenings in QCD-like theories. 
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In the AdS / CFT correspondence, the expectation value of the Wilson loop^ is related 
to the string partition function [26] 



{W{C)) 



(4.2) 



where we integrate over all the string worldsheets E with the boundary condition dH = C 
at the AdS boundary and S^GiJ^) corresponds to the Nambu-Goto action for the string 
worldsheet 



S 



NG 



J drda^ det {g^.udaX^dpx''). 



(4.3) 



Where g^j_jj is the bulk metric and l/2TTa' is the string tension; x^{t,(j) is the location 
of the string worldsheet in this {d + l)-dimensional spacetime. In the strong coupling 
limit, a' ^ 1 and we can perform a saddle point approximation, yielding 



{W{C)) 



-S'jvg(So) 



(4.4) 



where Eq represents the minimal area surface with the boundary condition SSq = C. 

4.1 Rectangular Wilson loop: Finite temperature expansion 

In this section, we will consider spatial rectangular Wilson loops specified by x = a;^ G 
[— |, |] , y = E f 5 f] and t = = for i > 2 (see figure 2(b) for a pictorial 
representation). And we will assume that L ^ I. The two dimensional worldsheet can 
be parameterized by the coordinates a'^ = (t, cr). We will pick cr = r, t = y as the 
worldsheet coordinates. In the limit L — )■ oo, the extremal surface is translationally 
invariant along the y-direction. Therefore, only coordinate x = x{r) has a nontrivial 
profile and the induced metric on the worldsheet is given by 



ds'L = r^dy^ + 



2 /2 I 

r X + 



r2 ( 1 - ^ 



dr' 



(4.5) 



and the action (4.3) becomes 



Sng 



dr 



\ 



r^x'^ + 



1 - 



(4.6) 



The minimal area surface Eq can be found by solving the equation of motion arising 
from the above action 



^In these gauge theories the Wilson loop operators are typically given by some generalization of 
equation (4.1). 
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where, is an integral of motion and r = represents the point of closest approach of 
the extremal surface.^ Each surface, just like the geodesic case, has two branches, joined 
smoothly at (r = rc,x = 0) and Tc can be determined using the boundary conditions: 

x(oo) = (4.8) 



That leads to 



"2 



/ _ rldr / rjj^ ^'^ 



r\ [I 



, , , . . (4.9) 

Again, analogous to the two-point function case, we can develop a systematic expansion, 

^_ 1 - T[\ + n]T[\{^ + nd)] /r^y^ 

2r,^ T[l + n]V + nd)] \rj ' ^ ' ' 

For large n the series goes as ~ ^{rH/f'c)"''^ and the series converges as long as th/tc < 1. 
Again it can be shown [20] that at any finite temperature > and hence the sum 
(4.10) is well-defined. Next we will calculate the on-shell action 



2L dr ( _ r 



d \ -1/2 
H 



The factor of 2 appears because of the two branches of the extremal surface. The on- 
shell Sng is divergent;'' however, it can be renormalized by introducing a UV-cutoff 
and then subtracting the boundary term rfeL/(7ra') 

SNG;ren = SnG (4-12) 



Now a series expansion for SNG;ren can be obtained 
Lr, 



SNG:ren — , 



i^r[| + n]r[i(-i + nd)] (m'^'^'' 



r(i)^ 4^ r[i + n]r[|(i + nrf)] 



(4.13) 



Again it can be shown that the series (4.13) converges for rujr^ < \. In principle, 
the rest of the calculation is straight forward. We have to solve equation (4.10) for Tc 
and then use that Tc in equation (4.13) to get SNG;ren- The expectation value of the 
renormalized Wilson loop is then given by. 



^It should be noted that rj, is an observable dependent quantity and hence Tc for the spatial rectan- 
gular Wilson loop is not the same as Vc for equal-time two-point function. 

^Note that the divergent piece is independent of both temperature and dimension. 
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4.2 Rectangular Wilson loop: Low temperature limit 

At low temperature (in this context, low temperature means Tl <^ 1), the calculation is 
similar to the two-point function case (see appendix A). The leading contributions come 
from the boundary and hence we can solve equation (4.10) for Vc order by order, leading 
to finite temperature corrections to the expectation value of the Wilson loop. Solving 
(4.10), at first order in {IthY, we obtain 



2v^7r3/2 



1 + 



2 2 2 vr 2 " 



Equation (4.13) then tells us 



(M'^ + 0(/r^) 



2d 



(4.15) 



S 



An' 



NG ■,ren 



«'r(|)4 



L 



where Pi is a numerical constant given by 



(4.16) 



2 2 27r 2 " 



ir(i)-+^r(f-i) 



r(f + !) 



(4.17) 



'SnG:i 



Expectation value of the renormalized Wilson loop is then given by {W{C)) = e 
Constant Vi is negative indicating a decrease in the expectation value as the system 
is heated. Replacing rn by corresponding temperature (2.5), for the d-dimensional 
boundary theory we obtain 



S 



An' 



NG 



a' 



\ fL\ \ /AnTlV fAnTl\ 

)(t) ) ) 



2d 



(4.18) 



where, we have restored AdS radius R by dimensional analysis. Particularly, for ci = 4 
at low temperature limit, we get 
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(4.19) 



4.3 Rectangular Wilson loop: High temperature limit 



At high temperature {thI ^ 1), just like the geodesic approaches th and the 

leading contribution comes from the near horizon region. We can rewrite equation (4.13) 
in a way so that we can take the limit — ?■ r// 



S 



NG:ren 



Lrr 



na 



2V2n^/' 



2 2^nd-l 



n=l 



r[l + n]r [\{5 + nd)] 



(4.20) 
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It is easy to check that the infinite series in the last equation converges for > and 
we can safely take the limit — )■ 



Lr 
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(4.21) 



2 ' 2^nrf-l T[l + n]T [1(5 + nd)] 

Therefore, at high temperature the expectation value of the Wilson loop is approximately 
given by 

{W{C)) 

where, A = IL is the area of the Wilson loop and 



(4.22) 
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In terms of temperature T, the asymptotic behavior of the expectation value of the 
Wilson loop in ci— dimensions is given by (after restoring AdS radius R) 



{W{C)) ^ exp 



a' 
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c/2 
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(4.24) 



At large temperature, the term proportional to the area in SNG;ren dominates indicating 
an area law which is consistent with [28]. The area term comes from the near horizon 
part of the extremal surface and could be understood easily by making an argument 
similar to the geodesic case (Figure (1) could be think of as a section of the extremal 
surface). The term independent of / receives contributions from the full bulk geometry. 

At high temperature {Tl ^ 1), similar to the two-point function case, approaches 
rn exponentially fast (see appendix B) 



(4.25) 



where S^i is a numerical constant given by, 
1 



d 



exp 



^^2V2ir^^lj^(T[^ + n]T[\{3 + nd)] 2 
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(4.26) 



Using this expression for rc, it has been shown in appendix B that the correction to the 
high temperature result (4.24) is very small 



_ f^'^\ fSTrAT^ ALTWd 4:LT 
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where the dots represent the higher order correction terms. 
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5 Entanglement entropy 



The entanglement entropy in quantum field theories or quantum many body systems is 
another important non-local quantity that can be employed to probe quantum properties 
of the system. Consider a quantum field theory with many degrees of freedom at zero 
temperature. Let us assume that the system is described by the pure ground state |\E'), 
which does not have any degeneracy. The density matrix of the state is 

Ptot = |^)(*| (5.1) 

and hence the von Neumann entropy of the total system, defined as Stot = — tr(ptot In Ptot) = 
0, does not contain any useful information. On the other hand, the entanglement en- 
tropy is non-vanishing even at zero temperature and several aspects of quantum many 
body physics can be understood by studying the entanglement entropy associated with 
volumes of different shapes and sizes. 




(a) (b) 

Figure 2. (a) The total system can be divided into two subsystems A and B; the entanglement 
entropy Sa measures the amount of information loss because of smearing out in region B. (b) 
A schematic diagram of the rectangular geometry and the corresponding extremal surface used 
for the calculations of spatial Wilson loops and the entanglement entropy. 



Now consider dividing the total system into two subsystems A and B and imagine 
an observer who has access only to the subsystem A (see figure 2(a)). The total Hilbert 
space is a direct product of two spaces Titot = 'Ha ®'Hb- The observer who is restricted 
to A will describe the total system by the reduced density matrix 

Pa = tiBptot, (5.2) 

where, degrees of freedom inside the subsystem B have been traced out. The entangle- 
ment entropy of the subsystem A is now defined as the von Neumann entropy of the 
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reduced density matrix pA 

Sa = -tr a{pa'^'o. Pa). (5.3) 

The entanglement entropy 5*^ describes the amount of information loss because of smear- 
ing out in region B and provides us with a convenient way to measure how the subsystems 
A and B are correlated with each other. At finite temperature, the entanglement entropy 
is defined in the same way but with the thermal density matrix p = e~^^'^, where H is 
the total Hamiltonian. 

A precise prescription for computing entanglement entropy for strongly coupled field 
theories with AdS duals was proposed in [30] and later generalized in [31] (for a good 
review see [32]). According to the proposal, the entanglement entropy Sa is given by 



Area (7^) 



4G 



(d+i) 

N 



(5.4) 



where, Cj^'*'^^ is the ((i + l)-dimensional Newton's constant, ja is the (d — l)-dimensional 
minimal area surface in the bulk whose boundary is given by the boundary of the region 
A: d'jA = dA. The area of the surface ■ja is denoted by Area (7a)- 

From this prescription it is very clear that for d = 2 and 3, the entanglement 
entropy calculation is the same as the equal time two-point function calculation and the 
spatial Wilson loop calculation respectively. Therefore, only for d > 4, the entanglement 
entropy contains non-trivial information. 

5.1 Entanglement entropy: Finite temperature expansion 

In this section we will compute the entanglement entropy for a strip (see figure 2(b) for 
a schematic diagram) specified by 
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(5.5) 



with L — 7- 00. Extremal surface is translationally invariant along x*, z = 2, ci — 2 and 
the profile of the surface in the bulk is x{r). Area of this surface is given by 



A = L 
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(5.6) 



This action leads to the equation of motion 

dx ^ 
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1 - 



(5.7) 



where, is an integral of motion and r = represents the point of closest approach 
of the extremal surface (again it should be noted that is an observable dependent 
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quantity). Similar to the previous cases, each surface has two branches, joined smoothly 
at (r = rc,x = 0) and can be determined using the boundary conditions: 



That leads to 
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(5.8) 
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By now it is obvious that we will do an expansion 
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Again it can be shown that the series converges for Vc > th- The area of the extremal 
surface is given by. 
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(5.11) 



This area is infinite indicating that the entanglement entropy has a divergence. In a 
field theory the entanglement entropy is always divergent because there are too many 
degrees of freedom. Therefore, we can write 



A — Afiiy + A 



div ~r ■'^finite- 



(5.12) 



The divergent piece is temperature independent and hence easy to compute. We will 
introduce an infrared cut off r^, which corresponds to the ultraviolet cut off a = 1/rf, (or 
a lattice spacing) of the boundary theory^ 
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(5.13) 



The divergence is proportional to the area of the boundary of A which is expected 
since the entanglement between A and B is strongest at the boundary dA. This area 
law behavior of the divergent piece is well understood from field theory computations 
[33-38]. 



'We are working with AdS radius R = I. Restoring R, the lattice spacing is given by a = 
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Analogous to the previous cases, we can do an expansion {d ^ 2) for the finite part 
of the area^ 
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(5.14) 



It can be shown that this series converges for > r^- Now, the rest of the procedure is 
simple and familiar. We have to solve equation (5.10) for and then we can calculate 
area by using equation (5.14). Entanglement entropy of the rectangular strip can be 
computed using the relation (5.4). In practice, this procedure can not be performed 
analytically at finite temperature. However, we can extract low and high temperature 
behavior of the entanglement entropy from equations (5.10,5.14). 

5.2 Entanglement entropy: Low temperature limit 

The temperature now should be measured with respect to ~ 1//; therefore, low temper- 
ature means Tl <^ 1. At low temperature, ^ and the leading contributions to the 
area come from the boundary which is still AdS. Therefore we should expect the zero 
temperature entanglement entropy as the leading term. Finite temperature corrections 
correspond to the deviation of the bulk geometry from pure AdS. At low temperature, 
the extremal surface is restricted to be near the boundary region and hence the devia- 
tion is small and can be computed perturbatively. At low temperature limit {thI <^ 1), 
equation (5.10) can be solved for and at first order in [vhIY, we obtain 
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(5.15) 



Now using equation (5.14), at first order in {thIY, we get (the calculation is similar to 
the two-point function calculation explained in appendix A) 
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^We will not consider the entanglement entropy for d = 2 case because the calculations are exactly the 
same as the two-point function case and hence has been studied analytically. In c? = 2, the entanglement 
entropy has a logarithmic divergence [39, 40], indicating a violation of the simple area law. 
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where, numerical constants Sq,Si are given by 
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Therefore, following equation (5.4), after restoring AdS radius i?, the entanglement en- 
tropy of the rectangular strip for the d-dimensional boundary theory at low temperature 
{Tl ^ 1) is given by. 
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(5.19) 

Note that in the limit T — )■ 0, we recover the well known results of [42]. Particularly for 
d = 4 we have, 
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(5.20) 



5.3 Entanglement entropy: High temperature limit 



At high temperature (i.e. Tl ^ 1), it is not very difficult to find out the asymptotic 
behavior of the entanglement entropy. However, it is more difficult to do a systematic 
expansion. At high temperature, the extremal surface tends to wrap a part of the horizon 
and the leading contribution comes from this near horizon part of the surface. For the 
subleading terms, the full bulk geometry contributes and they are more interesting. 

At very high temperature [thI ^ 1), Tc approaches th- By now we know how to 
determine the high temperature asymptotic behavior; we will rewrite equation (5.14) in 
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a way that allows us to take the limit — )■ th without encountering any divergence. 
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(5.21) 



The infinite series in the last equation for large n goes as ~ ^{th /tc)^'^ and thus the 
limit Tc — th exists. At high temperature ~ and the leading behavior can be 
determined by taking the limit — )■ th in the last equation 
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where, Shigh is another numerical constant given by 
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Hence, the entanglement entropy of the rectangular strip for the rf-dimensional boundary 
theory at high temperature is given by. 
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(5.24) 



where V = IL'^"'^ is the volume of the rectangular strip and R is the AdS radius. 

The divergent part of the entanglement entropy is temperature independent and 
thus it does not contain any new information. The leading finite piece in equation (5.24) 
is proportional to the volume of the rectangular strip and it is just the thermal entropy 
of the region A. The extrinsic nature of the leading term at high temperature can be 
understood very easily by looking at the extremal surface for rnl ^ 1. In this limit, the 
extremal surface tends to wrap a part of the horizon and the actual U-shaped surface 
can be approximated by a surface that consists of a; = —l/2,r = r^^x = 1/2 (one can 
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think of figure (1) as a section of the extremal surface). At high temperature hmit, the 
most dominant contribution to the area of the extremal surface comes from the near 
horizon part which can be guessed from this approximate surface. The area of the near 
horizon part of the approximate surface is A ~ r'jj . Therefore, it is expected that 
the leading term goes as Sa ~ VT'^~^. On the other hand, the other term ~ {LTY~'^ 
is more interesting. This term corresponds to the entanglement between the region A 
and the outside and it is proportional to the area of the boundary OA of A because 
the entanglement is strongest at the boundary. One can guess the functional form of 
this term from the approximate surface. However, the value of the numerical constant 
Shigh obtained from the approximate surface is inaccurate. Although the actual extremal 
surface approaches the approximate one as thI is increased, it can be shown that the 
area of this approximate surface is always greater than the actual extremal surface [41] 
and even in the limit thI — oo, two do not coincide. 

So far we have used the fact that at high temperature, ~ th- [20] shows that the 
extremal surface approaches the horizon but it always stays finite distance above the 
horizon. It is an interesting exercise to see exactly how fast approaches th. As shown 
in appendix (B), in the limit thI ^ 1 
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where £ent is a constant given by. 
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Using equation (5.25), we can calculate the next order correction to the high temperature 
entanglement entropy (for details see appendix (B)). The subleading term is exponen- 
tially suppressed 
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(5.27) 



^2d{d - 1) A T'^'^ exp |-v/(rf- l)/(2rf) 47rT/| + . 

where A is the area A = 2L'^^'^ and the dots represent the higher order correction terms. 

5.4 Entanglement entropy of a generic region at high temperature 

The calculation of the entanglement entropy of an infinite rectangular strip suggests that 
the general form of the finite part of the high temperature answer does not particularly 
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depend on the shape. One expects that the finite part of the entanglement entropy of a 
region A for a d— dimensional {d > 2) boundary theory with AdS-dual should be given 
by 



SA-jinite = Co [T"*"^ Volume(A) + ci T'^"^Area(M)] + sub-leading terms, (5.28) 

provided the temperature T ^ l/l, where / is the smallest length scale of the region A}^ 
Co is a constant that depends on the actual theory and Ci is a constant that depends on 
the shape of the region A. The first term represents the thermal entropy of the region 
A. The second term is proportional to the area of the boundary of A and corresponds 
to the entanglement between region A and the outside. 

6 J\f = 4: super- Yang-Mills in (3 + 1) dimensions at finite tem- 
perature 

In this section, we will use all the tools we have developed in the previous sections to 
study the prototype case of = 4 super- Yang-Mills with gauge group U{N) at finite 
temperature. The AdS / CFT correspondence relates this theory to type IIB string theory 
on asymptotically AdSs x spacetime. And as a consequence, in the limit ^ 1, A = 
9ym^ ^ 1) theory can be well approximated by the classical supergravity. 

Af = 4 super- Yang-Mills theory is the unique maximally supersymmetric gauge 
theory in (3-1-1) dimensions and it consists of gauge field A^, six scalar fields = 1, 6 
and four Weyl fermions Xa, all in the adjoint representation of U{N). The Euclidean 
action for the theory has the schematic form 



Now at finite temperature, in the Euclidean description, the system lives on x S*^, 
the circle direction corresponds to the Euclidean time with period T^^. At length scale 
/ ^ at least for weak coupling, one can perform a Kaluza-Kline reduction along 
the circle. All the fermions of the theory get a mass of order T at tree level because of 
the antiperiodic boundary conditions around the circle. The scalars get a mass of the 
order g^T at one loop level through their couplings to fermions. Therefore, at least for 
weak coupling, at large distance / 2> 1/T, the theory reduces to non-supersymmetric 
Yang-Mills theory in three dimensions. Because the theory has a mass gap, at large tem- 
perature (Tl 1), it is expected that correlation function {0{x,t)0{y,t)) ~ e"'^"^'"^. 

^''One can make a similar argument for the spatial Wilson loop. 



s = - 




(6.1) 
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Things are more complicated at strong coupling. Fortunately, thanks to AdS/CFT, 
at strong coupling, we can use the supergravity solution to study the system. The mass 
of the excitations of the QCD strings [28] is Mg = X^^^T and the mass associated with 
the compactification is Mc = T . In the strong coupling region (A ^ 1), Mg ^ and as 
a result QCD strings can "see" the compactified circle. Therefore at high temperature 
{Tl ^ 1), the theory does not reduce to pure Yang-Mills in three dimensions but = 4 
super- Yang-Mills in four dimensions compactified on a circle. 

Let us now put h back in the picture for clarity. The action (6.1) has a dimension 
ML. The equal time two-point function of the gauge invariant variable O = ^TrF"^ is 
a useful non-local observable to understand the low and high temperature behavior of 
the theory. The operator O has conformal dimension A = 4 with a unit ML~^}^ Now 
using equation (3.28) and the relation 
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at low temperature {T\x — y\ h) we obtain 
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The factor of h"^ is there because {0{x,t)0{y,t)) has dimension M'^L~^. The factor of 
comes from comparing the result with standard zero temperature result. Similarly, 
at high temperature (T\x — y\ ^ h), using equation (3.37), we obtain 

{0{x,t)0{y,t)) = ^<^^lf'^\ -^-T\.-y\/n. ^ 2v^^4 g-^-^l^-^l/''' + ...1 . (6.4) 

Where = 4.51 and £4 = 2.75. In the last equation, again we have used dimen- 
sional analysis to recover h. Indeed, at high temperature, the two-point function decays 
exponentially. 

Expectation values of the Wilson loops are another set of gauge-invariant non-local 
observables that are useful to understand the non-perturbative behavior of non-Abelian 
gauge theories. In this Af = 4 SYM theory, one should generalize the Wilson loop 
expression (4.1) in the following way 

W{C) = ^Tr (pe^cH^^-'+^-f^)^ , (6.5) 

where, x^{s) parametrizes the path C and n is a unit vector in the {(pi} space. Expecta- 
tion value of a rectangular infinite Wilson loop with long side of the loop extends along 



"'^"'^The saddle point approximation that we have used to obtain the two-point function is a good 
approximation only for large A. For definiteness we will use O = ^TrF^ just as an example and 
consider only the geodesic contributions and ignore the sub-leading corrections. 
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Figure 3. Variation of f{x,y) = ^ln{OA{t,x)OA{t,y)) with T\x — y\ for the 4-dimensional 
A/" = 4 SYM theory. The sohd black Une represents the exact numerical result. Blue and red 
lines represent the two-point functions computed using low and high temperature approxima- 
tions, respectively. 



the time direction gives potential energy between a static quark-antiquark pair. Unlike 
the spatial Wilson loop case a sharp transition takes place at T ~ 1// for this case [27]. 
At low temperature, the U-shaped solutions exist but at high temperature no nontriv- 
ial solutions exist. At high temperature, the extremal surface consists of two disjoint 
vertical surfaces ending at the horizon. Therefore, at high temperature the Wilson loop 
is independent of separation I. In contrast, the transition is more gradual in the case 
of spatial rectangular (infinite) Wilson loops. In A/" = 4 SYM theory, the expectation 
value of a spatial rectangular (infinite) Wilson loop is given by. 
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and using equation (4.18) at low temperature {Tl ^ h), we obtain 
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At high temperature {Tl ^ h), equation (4.27) leads to 
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(6.8) 



where, A = IL is the area of the loop and E^i = 1.66. We have used the fact that for 
this theory /a' = a/A. 
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Figure 4. Variation of SNG;ren (in the units of -^^5^) of the spatial rectangular (infinite) 
Wilson loop with Tl for the 4-dimensional A/" = 4 SYM theory. The expectation value of 
the Wilson loop is given by {W{C)) = e~^^'^''^^^. The solid black line represents the exact 
numerical result. Blue and red lines represent the two-point functions computed using low and 
high temperature approximations, respectively. 



In a field theory the entanglement entropy is always divergent because there are too 
many degrees of freedom and we can write 



finite- 



(6.9) 



We will consider the entanglement entropy of an infinite rectangular strip (see section 
5). For the case in hand the divergent piece is temperature independent 
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Where we have used the AdS/CFT dictionary to write R^/4G 
temperature (Tl <^ h), the finite part is given by (see equation (5.19)) 
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And at high temperature [Tl ^ h), following equation (5.27) we obtain 
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where, £f.nt = 1-174, V = IL^ is the volume of the rectangular strip and A = 2L^ is the 
area. Note that the leading behavior is consistent with the observation of [42]. 



- 26 - 



'^finite 




Figure 5. Variation of SfinUe^ the finite part of the entanglement entropy (in the units of 
) of an infinite rectangular strip with Tl for the 4-dimensional Af = 4 SYM theory. The 
solid black line represents the exact numerical result. Blue and red lines represent the two-point 
functions computed using low and high temperature approximations, respectively. 

7 Conclusions 

In this article, we have studied the high and low temperature behavior of some non-local 
observables in strongly coupled large-N gauge theories. At low temperature, the leading 
term is obviously the zero temperature contribution and the sub-leading terms indicate 
a decrease in correlations as the system is heated. At high temperature, the leading 
contributions come from the near horizon part of the bulk. The full bulk geometry con- 
tributes significantly to the sub-leading terms. We have investigated the contributions 
from different regions of the bulk for the equal time two-point function, rectangular (in- 
finite) spatial Wilson loop and entanglement entropy of an infinite rectangular strip. In 
the case of spatial Wilson loops and the entanglement entropy, our calculations suggest 
that the general form of the high temperature answer does not particularly depend on 
the shape. In particular, it is interesting that at high temperature we observe an area law 
behavior of the finite subleading term of the entanglement entropy. This is consistent 
with what is expected from Quantum Information Theory when the two-point function 
decays [43]. However, the Quantum Information Theory argument is qualitative and 
can not be extended easily to non-zero temperature. 

In forthcoming work [44] we are generalizing this approach to include finite densities. 
We also plan to extend our techniques to the analysis of holographic thermalization which 
has attracted much attention in recent years [45-50]. An analytic understanding of ther- 
malization at different limits can possibly provide new insight into the non-equilibrium 
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physics. 
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A Two-point function: Low temperature expansion 

At low temperature, ru/rc ^ 1. From equation (3.13), keeping only few subleading 
terms, we get 
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Solving this perturbatively we obtain 
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In this low temperature limit, equation (3.18) becomes 
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Now using equation (A.2), we get 
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Therefore, the two-point function is 
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where again / 
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B High temperature subleading term 
B.l Two-point function 

We have mentioned that at high temperature Tc ^ th- Now we will try to figure out how 
fast Vc approaches rn- For large n, the series (3.13) goes as ~ ^ {rn/rc)"''^ and hence it 
diverges for Vc = th- We can isolate the divergent piece from the infinite series (3.13) , 
yielding 
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The infinite series in the last equation now converges even for = th- It is convenient 
to write = r//(l + e) for the purpose of solving the last equation. It is very clear 
from previous discussions that at high temperature e <^ 1. Therefore, from the last 
expression, we obtain 
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Solving the last equation for e, at the leading order we get 
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Now we will compute the next order correction to the high temperature result (3.32). 
First we will start with equation (3.30) 
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We can not yet write = r//(l + e) and do an expansion for small e. The term linear 
in e contains an infinite series that does not converge. Therefore, in order to obtain a 
finite first order expression in e, we will write the last equation in the following way 
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where, Li is the PolyLog function. We have summed the part of the infinite series that 
could diverge in the first oder in e. Now, we can write Tc = r/^(l+e) and do an expansion. 
In the first order in e we obtain, 
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In the first line, we have all the terms that contribute in the leading order at high 
temperature (we will call it Cieading) and in the second line we have all the subleading 
terms. We can simplify the second line farther by using equation (B.3) 



^ren ^leading 

Therefore, using (3.4) we finally obtain 
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where Sd is given by equation (B.4) and the dots represent the higher order correction 
terms. 



B.2 Rectangular Wilson loop 

Again we will isolate the divergent part of the series (4.10) 
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It is again convenient to write Tc = rni^ + e) and at high temperature e ^ 1. Therefore, 
from the last expression, we obtain 
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Solving the last equation for e, at the leading order we get 
where Eu,i is a constant given by, 
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Now we will compute the next order correction to high temperature result (4.24). 
First we will start with 
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Before, we write Tc = rH{^ + e) and do an expansion for small e, in order to obtain a 
finite first order expression in e, we will write the last equation in the following way 
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Where, Li is the PolyLog function. We have again summed the part of the infinite series 
that could diverge at first order in e. Now, we will write Vc = r//(l + e) and do an 
expansion. At first order in e we obtain. 



S 



Lr 



H 



NG:ren 



Tra' 



+ 



eLr 



H 



' T{\Y + 2 +2^nrf-l V[l + n]V[\{h + nd)\ 
2y27r3/2 1 



iy- /^r[i+n]r[i(3 + nrf)] 2 



n=l 



r[l + n]r [\{b + nd)\ 



d^/^n 



(B.15) 



- 31 - 



In the first line, we have all the terms that contribute in the leading order at high 
temperature (we will call it Sjyc-Mading)- We can simplify the subleading terms farther 
by using equation (B.ll) 
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Therefore, finally we have 
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B.3 Entanglement entropy 



We know that at high temperature ~ th. Now again we will try to figure out how 
fast Tc approaches th- For large n, the series (5.10) goes as ~ ^ {jhI'^cT'^ and hence it 
diverges for = th- We can isolate the divergent piece from the infinite series (5.10), 
obtaining 



1 




r[i + n]r 


2(d-l) 


r[i + n]r 


dn+1 
2{d-l) 




In 




l)d 



2^{d- l)d n 




rH_ 



(B.18) 



The infinite series in the last equation now converges even for Vc = th- Analogous to the 
last two cases, now it is convenient to write Tc = r/^ (1 + e). It is very clear from previous 
discussions that at high temperature e -C 1. Therefore, from the last expression, we 
obtain 
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Solving the last equation for e, at the leading order we get 
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where 8ent is a constant given by, 
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Now we will compute the next order correction to the high temperature result (5.22). 
First we will start with equation (5.21) 
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We will write = r//(l + e) and do an expansion for small e and then we will use 
equation (B.20). But before that we should note that the term linear in e contains an 
infinite series that does not converge. We will isolate and resum part of the last equation 
that could diverge at first order in e, in the following way 
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Where, Li is the PolyLog function. Now, we will write Tc = ri^(l + e) and do an 
expansion. In the first order in e we obtain, 

1 



A 



finite 



1 + 



Shi 



high 



+ 2L'^-V^-2e 



-(rf-i)E 

n=l 



7F(rf- i)r 



2[d-l) 



r 



1 



2{d-l) 



lrH + 



1 



T[\+n]V 


' d{n+l)' 
2{d-l) 


T[l+n]T 


dn+1 
2{d-l) 



2d 
1 



-l + log(ed)) 



V2^{d-l)d 



n 



(B.24) 



- 33 - 



In the first line, we have all the terms that contribute in the leading order at high 
temperature and Shigh is given by equation (5.23). We can simplify the subleading 
terms farther by using equation (B.20) 
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Therefore, finally we have 
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